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Abstract 

In this letter we discuss the Non-gaussian statistics considering two aspects. In the first, we 
show that the Maxwell's first derivation of the stationary distribution function for a dilute gas can 
be extended in the context of Kaniadakis statistics. The second one, by investigating the stellar 
system, we study the Kaniadakis analytical relation between the entropic parameter re and stellar 
polytrope index n. We compare also the Kaniadakis relation n = n(re) with n = n(q) proposed in 
the Tsallis framework. 
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I. INTRODUCTION 



The nonextensive statistical mechanics (NSM) [l| and extensive generalized power-law 
statistics jsj are based on the mathematical generalization of the Boltzmann exponential 
distribution, namely: 

/ ~ exp(— Energy /Thermal energy). (1) 



Whereas in the power-law distribution, / ~ [1 — (1— g) Energy /Thermal energy] 1 ^ 1 "^, 
obeys the Tsallis distribution of NSM jl| , the so-called ^-statistics 3 provides a power-law 
distribution function and the K-entropy which emerges in the context of the special relativity 
and in the kinetic interaction principle. Formally, the K-framework is defined by considering 
the expressions 

exp K (/) = (v / l + ^ 2 / 2 + ^) 1/K , (2) 

r - f~ K 



2k 



where the K-entropy associated with the K-statistics is given by 



S K 



d 3 pfln K f = -(ln K (f)). 



(3) 



(4) 



The expressions above reduces to the standard results in the limit k = 0. 

The ^-statistics leads to a class of one parameter deformed structures with interesting 
mathematical properties 3 , as well as a connection with the generalized Smoluchowski equa- 
tion yf and the relativistic nuclear equation of state for nuclear matter 5fl. A fundamental 

n 

test, the so-called Lesche stability was also studied in the context of K-framework [6j]. It 
was also shown that it is possible to obtain a consistent form for the entropy which is linked 



with a two-parameter deformations of 
and Kaniadakis logarithm behaviours 



ogarithm function and generalizes the Tsallis, Abe 
7]. There exist some applications related with the 

3| , quark-gluon plasma 



^-statistic, namely, cosmic rays flux, rain events in meteorology 
kinetic models describing a gas of interacting atoms and photons [9j, fracture propagation 
phenomena 10[, and construct financial models ll|. In the theoretical front, some studies 



on the canonical quantization of a classical system has also been investigated 
as if -theorem in relativistic and non- relativistic domain 13] 



121 ]. as well 



In the astrophysical domain, the first application has been the simulation in relativistic 
plasmas. In this regard, the power-law energy distribution provides a strong argument in 



2 



favour of the Kaniadakis statistics [14] . Additionally, the viability of the non-gaussian statis- 



tics has been investigated in stellar astrophysics viewpoint: the distri 
rotational velocity measurements of stars in the Pleiades open cluster 



jutions of projected 



sequence field stars 



15] and in the main 



16], as well as the strong dependence between the stellar-cluster ages 



and the power-law distributions 17]. 

The aim of this letter is twofold. First, to show that the Maxwell's first derivation 
of the stationary distribution function for a dilute gas can be extended in the context of 
Kaniadakis statistics; Second, considering the principle of maximum entropy for a stellar 
self-gravitating system, we investigate an analytical relation between entropic parameter k 
and stellar polytrope index n. It is also show that the function n = ti(k) has a similar 
behaviour to Tsallis expression n = n(q) 18]. 

This letter is organized as follows. In Section 2, we show the correspondence between the 
^-statistics introduced by Kaniadakis and the velocity distribution for a Maxwellian gas, by 
assuming a Non-gaussian generalization of the separability hyphothesis originally proposed 
by Maxwell 19]. In Section 3, we discuss a connection between the Kaniadakis statistics and 
polytropic index in the context of the self-gravitating system and we compare our results 



with ones studied in Ref. 



18] . Finally, Section 4 is devoted to conclusions and discussion. 



II. NON-GAUSSIAN MAXWELLIAN DISTRIBUTION FUNCTION 



In this section, in order to introduce the generalization of the Maxwell distribution in 
the context of Kaniadakis framework, let us consider a spatially homogeneous gas, supposed 
in equilibrium (or in stationary state) at temperature T, in such a way that F(v)d 3 v is 
the number of particles with velocity in the volume element d 3 v around v. In Maxwell's 
derivation, the 3-dimensional distribution is factorized (lottery assumption) and depends 



only on the magnitude of the velocity 



19 



20] 



F(v)d 3 v = f(v x )f(v y )f(v z )dv x dv y dv z 



(5) 



where v = ^Jv* + v % + v 2 z and F(v) is the standard Maxwellian distribution function, given 
by 

F(v) = Aexp(-mv 2 /k B T), (6) 

\ 3/2 



where A 



\2-Kk B T ) 



is the normalization constant. 



As a matter of fact, in the ^-statistic context described by (J4]), the starting basic hypoth- 
esis ([5]), which takes into account the isotropy of all velocity directions, must somewhat be 
modified. In statistical viewpoint, Maxwell's ansatz assume that the three components of the 
velocity are uncorrelated. However, this property does not hold in the systems endowed with 
long range interactions, or statistically correlated, where Kaniadakis distribution has been 



observed 



15 



161 ] . Notice that the Maxwell ansatz is equivalent to express InF as the sum 



of the logarithms of the one dimensional distribution functions associated with each velocity 
component. A simple and natural way to generalize this procedure within the Kaniadakis 
framework, in order to introduce correlations between the velocity components, would be to 
replace the logarithm function by a power law. However, in order to recover the ordinary 
logarithmic ansatz as a particular limiting case, it is convenient to express the power gener- 
alization in terms of the function ln re defined by Eq. ([3]), which is a combination of a power 
function plus appropriate constants. Mathematically, the consistent ^-generalization of ([5]) 
is given by 

F [yjvl + vl + vij d 3 v = exp K [ln K f(u x ) + hx K f(v y ) + 

ln K f(v z )]dv x dv y dv z , (7) 

where the K-exponencial and K-logarithm are given by identities (T5]) and ([3]). In particular, 
in the limit k = the standard expression ([5]) is recovered. Note also that ln K [exp tt (/)] = 
exp K [ln K (/)] = /, and dXn ^ = — y — ^ are satisfied. The logarithmic derivative of ([7]) with 
respect to Vi is 

^\R K F^Jvl + vl + vl) = -^ki K {exp K []n K f(v x ) + 

ln K f(v y ) + ln K f(v z )}}, (8) 

with i = x,y,z. 

Using the properties above mentioned we can write 

F K + p-K Q ( r Q 



2F dvi 

Equivalent ly, 



F (y/v* + v* + v*) = ^K/K)] • (9) 



7TB = -^~ ln «/ ( v i) . 10 

2F x ovi 



where x — \J v x + v y + v l- Now, defining 



$(x) = [X) , (11) 

yA) 2F x 
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the equation ({TO]) can be written as 



*(x) = --^K/(«01 

Vi dvi 



(12) 



The second member of Eq. ( TT2]) has dependence only on Vi, with i = x, y, z. Hence, the 
Eq. (fTTj) can be satisfied only if all its members are equal to one and the same constant, not 
depending on any of the velocity components. Thus we can make $(x) = -m7, obtaining 

1 d 



where 7 = . 



Vi dv. 



[\n K f(vi)] = -7717 



(13) 



In general, we have 



which provides 



rwyvf 



where 1/Z is the ^-normalization constant. 

In order to calculate the so-called ^-normalization, let us introduce the expression 

where n (=1,2,3) is the number from degree of freedom. Here, considering n = 1, a 



(14) 



(15) 



and x = av 2 , we obtain 



1 r°° 

z = 4= I 



x exp K (— x)dx. 



(16) 



2k B T 



(17) 



Using the generalized gamma functions |2] 



x r exp K (—x)dx 



[1 + (r - 2\K\)\2K\ r ] 
[1 - (r - 1)\k\} 2 - k 2 



r ^ 

|2«| 



1 5 



(18) 



and after some algebra, we obtain 



1 
4 



nk B T \k\^/ 2 r (l^f 

[i-iM]r(^ + i 



(19) 
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It is easy to see that the standard Maxwellian result (1/Z) = (m/2iikBT) 1 / 2 is recovered for 
gaussian limit k = 0. 

In order to calculate the complete distribution, we insert the expressions ( !T5|) into of Eq. 
([7]) for obtain 

,2 



F<«) = ~ z 




Kmv 



(20) 



2k B T J 2k B T 

Using the expression (Tl6l) for n = 3, we can show that the normalization for the complete 
distribution is given by 



1 


{ m | ac \ 3 ^ 2 




s, r 






1 H 


— k 


z ~ 






2 . 



r J- + 2 

1 l \2k\ ^ 4 



1 

2k| 



(21) 



As expected, the ^-distribution above is isotropic meaning that all velocity directions are 
also equivalent in this generalized context. Here, we emphasize that in the limit k = both 
the normalization and the distribution function reproduces the standard Maxwellian ([6]). 

III. NON-GAUSSIAN FRAMEWORK AND STELLAR POLYTROPES 



As mentioned in the introduction, we shall discuss a connection between the Kaniadakis 
statistics and polytropic index in the context of the self-gravitanting system. Following a 
procedure considered in Ref. 18|, let us start with the Kaniadakis generalized entropy of 
index n of the form 



f 



l + K 



/ 



1-K 



2k 



-dtt, 



(22) 



where dQ = d 3 rd 3 v, and the parameter k ^ provides a possible generalization of the 
Boltzmann Gibbs entropy. The extremum entropy state can be derived by varying S K with 
respect to /. Using the Lagrange multipliers a and (3, the extremum solution subject to 
constraints M = J fdQ and E = K + U = \ J v 2 fdVl + \$ 4>dVt 1 is obtained from 



8S K - a5M - (35E = 



(23) 



1 The mass and energy total of the self gravitating system governed by the Vlasov and Poisson equations. 
For more details, see Ref. 181 
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which leads to 



2k 



[(1 + «)/»- (l _«)/-"]- « 



-'(H) 



= 0. 



(24) 



We have used the relation J 5<f)fdQ = J <f)5fdQ for derivation of above expression. Here, 
<p denotes gravitational potential, and the equation (123]) must be satisfied independently of 
the choice of 5f. Thus, we obtain the following distribution function 



f(r,v)=B 



1/k 



(25) 



where the constants B and 0o are defined by 



B 



k(3 



1 + K 



1/k 



00 (C) 



■>/ k 2 (— a — /3e) 2 + 1 + — aft 



/3k 



(26) 



Now, comparing the distribution 



f l25l) with the polytropic spheres distribution / ~ £ n ~3/ 2 
(e is the relative energy of a star) [2lj], we obtain the relation between the stellar polytropic 
index and the entropic index, given by 



n 



3 1 
2 + «' 



(27) 



In order to compare the expression n = with the Tsallis relation n = n(q), we 

present the expression discussed in Ref. [18] 



n 



3 1 
2 + g- 1 ; 



(28) 



which include the isothermal situation for n — > oo, i.e., the Gaussian limit q — 1. In order 
to guarantee the conservation of mass and energy, the Tsallis' parameter g should be larger 
than 9/7. 

In Figure 1 we display the dependence between the entropic indexes (k; q) and the poly- 
tropic index (n), given by Eqs. ( |27|) and ( |28|) . respectively. We show that the polytropic 
index diverge for the non-gaussian parameters k — 0, q — 1, i.e., the isothermal spheres 
separates the polytropes into two branches. For k > and g > 1, n ranges from oo to 3/2. 



From / ~ e 



n-3/2 



we see that e is decreasing function of the energy and, when n = 3/2, 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Parameters (q, k ) 



FIG. 1: Polytrope index n as a function of the Kaniadakis (k) and Tsallis (q) parameters. The vertical 



dotted lines mark the limits of validity of ^-statistical (— 1 < k < 



1) fl. 



the distribution function becomes a constant independent of energy. On the other hand, it 
is well known that for any astrophysical system, n should be positive and higher than 1/2 



2l|. Thus the intervals < q < 1 and — 1 < k < represent forbidden regions, since the 



values of n index tend to be smaller than 1/2 and negative. 

It is also worth observing that for k < — 1 and q < also provides n > 1/2, however 
this limits on k; q violates the validity of ^-statistics [2J and the non-extensive if-theorem 



23(|, respectively. As is well known, for the polytropic index n > 5, the density falls off so 



slowly at large radii that the mass is infinite 2l|. Figure 1 shows that only the branches 
of the right-side are physically significant for n < 5. Therefore, the physical values of the 
Tsallis parameter q should be larger than 9/7 2 , and the Kaniadakis parameter k should be 
constrained to interval of validity k£ [2/7; 1] represented by the lower dashed rectangle. 



2 Investigating a collisionlcss stellar gas in the context from Non-extensive Kinetic Theory, the bound 
q > 9/7 also was calculated in Ref. [22[ 
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IV. CONCLUSIONS 



In this work we have studied the Kaniadakis statistic based on the generalized Maxwellian 
formulation for the ^statistic and, as one application, we investigate the physical effect on 
the stellar polytropic system. In the first part, we conclude that there is a Kanidakis 
velocity distribution given by Eq. (20), that is uniquely determined by the requirements 
of (i) isotropy and (ii) a suitable generalization of the factorizability condition. From the 
physical viewpoint, this generalization introduces statistical correlations between the velocity 
components, when we replace the logarithm function by a power law (similar argument 



was considered in Refs. 



23 



24j for the non-gaussian generalization of molecular chaos 
hypothesis). In particular, Maxwell expressions are recovered in the gaussian limit, n — > 0. 

From an application perspective, we have shown that the expressions for Tsallis and 
Kaniadakis, given by Eqs. (1281) and (1271) . present similar behaviours. However, the astro- 
physical limit on the polytropic index 1/2 < n < 5, provides the constrains k e [2/7; 1] and 
q > 9/7 for the Kaniadakis and Tsallis parameters, respectively. It is worth mentioning that 
the gaussian limit k — 0, equivalent the Tsallis expression Eq. 02 8p for q — 1, reproduces 
Maxwellian isothermal spheres or, equivalently, n = oo. 
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